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Abstract 

The Kuramoto-Daido model, which describes synchronization phenomena, is a system 
of ordinary differential equations on iV-torus defined as coupled harmonic oscillators, whose 
natural frequencies are drawn from some distribution function. In this paper, the continuous 
model for the Kuramoto-Daido model is introduced and the linear stability of its trivial solution 
(incoherent solution) is investigated. Kuramoto's transition point K c , at which the incoherent 
solution changes the stability, is derived for an arbitrary distribution function for natural 
frequencies. It is proved that if the coupling strength K is smaller than K c , the incoherent 
solution is asymptotically stable, while if K is larger than K c , it is unstable. 



§ 1. Introduction 

Collective synchronization phenomena are observed in a variety of areas such as 
chemical reactions, engineering circuits and biological populations [16]. In order to 
investigate such a phenomenon, Kuramoto |9j proposed a system of ordinary differential 
equations 

N 

(1.1) ^ = Ui + ^ sin ^' ~0i),i = h---,N, 

i=i 

where 9i G [0, 2tv) denotes the phase of an i-th oscillator on a circle, w ; G R denotes 
its natural frequency, K > is the coupling strength, and where N is the number 
of oscillators. Eq.( ll.ll) is derived by means of the averaging method from coupled 
dynamical systems having limit cycles, and now it is called the Kuramoto model. 

It is obvious that when K = 0, 9i(t) and 0j(t) rotate on a circle at different 
velocities unless u>i is equal to Uj, and it is true for sufficiently small K > 0. On the 
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Figure 1. The order parameter for the Kuramoto model. 



other hand, if K is sufficiently large, it is numerically observed that some of oscillators 
or all of them tend to rotate at the same velocity on average, which is called the 
synchronization [T§1 [33] • If N is small, such a transition from de-synchronization 
to synchronization may be well revealed by means of the bifurcation theory [31 EED, d2] • 
However, if N is large, it is difficult to investigate the transition from the view point of 
the bifurcation theory and it is still far from understood. 

In order to evaluate whether synchronization occurs or not, Kuramoto introduced 



which gives the centroid of oscillators, where r, ip G R. It seems that if synchronous 
state is formed, r(t) takes a positive number, while if de-synchronization is stable, r(t) is 
zero on time average (see FigJT]). Based on this observation and some formal calculation, 
Kuramoto conjectured a bifurcation diagram of r(t) as follows: 

Kuramoto's conjecture 

Suppose that N — > oo and natural frequencies u>iS are distributed according to 
a probability density function g(u). If g(u) is an even and unimodal function, then 
the bifurcation diagram of r(t) is given as Fig[2] (a); that is, if the coupling strength 
K is smaller than K c : = 2/(^(0)), then r(t) = is asymptotically stable. On the 
other hand, if K is larger than K c , there exists a positive constant r c such that r{t) = 
r c is asymptotically stable. Near the transition point K c , the scaling law of r c is of 



Now the value K c = 2/(71(7(0)) is called the Kuramoto's transition point. See [10] 
and [E] for the Kuramoto's discussion. 

Significant papers of Strogatz et al. [T9l EDI [IS] partially confirmed the Kuramoto's 
conjecture. Though their arguments are not rigorous from a mathematical view point, 




(1.2) 




0((K — Kc) 1 / 2 ). 
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Figure 2. Typical bifurcation diagrams of the order parameter for the cases that (a) 
g(u) is even and unimodal (b) g(u) is even and bimodal. Solid lines denote stable 
solutions and dotted lines denote unstable solutions. 

almost all of them are justified as will be done in this paper. In [2D], they introduced 
the continuous model for the Kuramoto model and investigated the linear stability of a 
trivial solution called the incoherent solution, which corresponds to the de-synchronous 
state r = 0. They derived the Kuramoto's transition point K c = 2/(^(0)) and showed 
that if K > K Cl the incoherent solution is unstable in the linear level (i.e. nonlinear 
terms are neglected). When K < K Cl the linear operator T, which defines the linearized 
equation of the continuous model around the incoherent solution, has no eigenvalues. 
However, in [19], they found that an analytic continuation of the resolvent (A— T) _1 may 
have poles {resonance poles) on the left half plane, and they remarked a possibility that 
resonance poles induce exponential decay of the order parameter. In [15], the stability 
of the partially locked state, which corresponds to a solution with positive constant 
r = r c , is investigated in the linear level. 

Despite the active interest in the case that the distribution function g(u) is even 
and unimodal, bifurcation diagrams of r for g(u) other than the even and unimodal 
cases are not revealed well. Martens et al. [13] investigated the bifurcation diagram for 
a bimodal g(u) which consists of two Lorentzian distributions. In particular, they found 
that stable partially locked states can coexist with stable incoherent solutions if K is 
slightly smaller than K c (see FigJ2] (b)). Such a diagram seems to be common for any 
bimodal distributions. 

A simple extension of the Kuramoto model defined to be 



is called the Kuramoto-Daido model [H El El E], where the 27r-periodic function / : 
R — > R is called the coupling function. Daido [7] investigated bifurcation diagrams of 
the order parameter for the Kuramoto-Daido model with even and unimodal g{oS) by 



(1.3) 




3=1 
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a similar argument to Kuramoto's one. He found that if f(9) ^ sin (9, partially locked 
states may coexist with stable incoherent solutions even if g(uj) is even and unimodal. 

All such studies by physicist are based on formal calculations and numerical simu- 
lations. The purpose of this paper is to justify and extend their results as mathematics 
for the Kuramoto-Daido model with any distribution function g{uj). The continuous 
model for the Kuramoto-Daido model is introduced and the linear stability of the in- 
coherent solution is studied. In particular, the spectrum and the semigroup of a linear 
operator T, which is obtained by linearizing the continuous model around the incoherent 
solution, will be investigated in detail. At first, a formula for obtaining the transition 
point K c for an arbitrary distribution g(uj) is derived. As a corollary, the Kuramoto's 
transition point K c = 2/(ng(0)) is obtained if g(u:) is an even and unimodal function. If 
K > K c , it is proved that the incoherent solution is unstable because the operator T has 
eigenvalues on the right half plane. It means that if the coupling strength K is large, 
the de-synchronous state is unstable and thus synchronization may occur. On the other 
hand, if < K < K Cl it will be shown that the spectrum of the operator T consists 
of the continuous spectrum and it lies on the imaginary axis. Thus the stability of the 
incoherent solution is nontrivial. Despite this fact, under appropriate assumptions for 
g{u>), the order parameter proves to decay exponentially because of existence of reso- 
nance poles on the left half plane as was expected by Strogatz et al. [TH]. It suggests 
that in general, linear stability of a trivial solution of a linear equation on an infinite 
dimensional space is determined by not only the spectrum of the linear operator but 
also its resonance poles. 

§ 2. Continuous model 

In this section, we introduce a continuous model of the Kuramoto-Daido model and 
show a few properties of it. 

Let us consider the Kuramoto-Daido model (11.31) . We suppose that the coupling 
function / is a C 1 periodic function with the period 2-k. It is expanded in a Fourier 
series as 

oo 

(2.1) /(*)= f L eC. 

l= — oo 

We can suppose that fo — without loss of generality because fo is renormalized into 
the constants Ui. For the Kuramoto model (/(#) = sin6>), f±i = ±l/(2\/— 1) and 
fl = 0{l^ ±1)- Following Daido [7], we introduce the generalized order parameters Z® 
by 

1 - 

(2.2) z° k (t) = -^Y, e jit) > * = 0,±1,±2,.... 
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In particular, Z\ is the order parameter defined in Section 1. By using them, Eq. (jl.3j) 
is rewritten as 



(2.3) 



dfk 
dt 



l=— oo 



Motivated by these equations, we introduce a continuous model of the Kuramoto-Daido 
model, which is an evolution equation of a probability density function p t = pt(9, u) on 
S 1 x R parameterized by t G R, as 



(2.4) 



dt 



a I oo 



-iw 



z° k (t) 

RJO 

k Po(0,u) = h(9,u) 




2n 



-ike 



Pt(9, uj)d6d<jj, 



where h(6,u) is an initial density function. The Z®(t) is a continuous version of Z^(t), 
and we also call it the generalized order parameter. We can prove that Eq. (12.41) is proper 
in the sense that Zf&t) -> Zf&t) as — > oo under some assumptions, although the proof 
is not given in this paper. If we regard 



v := u 



+ K J2 fiZ?(t)e 



-ue 



l= — oo 



as a velocity field, Eq. (|2.4p provides an equation of continuity dpt/dt + d(ptv)/d9 = 
known in fluid dynamics. It is easy to prove the low of conservation of mass: 



(2.5) 




2tt 



Pt(0,u)d9 




2tt 



h(6,u)d6=: g{u) 



A function g defined as above gives a probability density function for natural frequencies 
iv G R such that J R g(u)du = 1. 

By using the characteristic curve method, Eq. (j2.4p is formally integrated as follows: 
Consider the equation 



(2.6) 



dx 



l= — oo 



which defines a characteristic curve. Let x = x(t,s;9,uj) be a solution of Eq. (j2.6|) 
satisfying x(s, s; 6, u) = 9. Then, p t is given as 



(2.7) p t (9,uj) = h(x(0,t;d,u),u)exp K / ^ V=Ti/ I Z,°(s)e 



t 



-iix{s,t-e,uj) 



ds 
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By using Eq. fl2. Tj) , it is easy to show the equality 

/» /»27T i> /»27T 

(2.8) / / a(6,u)p t (6,uj)d6duj = / a(x(t,0;d,u),u)h(d,u)<ffldu), 



jujo JnJo 

for any continuous function a(9,u). In particular, the generalized order parameters 
Zfc(t) are rewritten as 

(2.9) Z° k (t)= [ [ e V=Ikx ^°' e ^h(9,u)d9(ko. 

JnJo 

Substituting it into Eqs. (12.61) . ( 12 .7p . we obtain 

(2.10) ^-x(t,s;9,u) = u + K I [ f{x{t,Q\6' ,uj') - x{t,s\B,u))h{e' ,u}')d9'duj' , 
dt JnJo 

and 

p t (9,uj) = h(x(0, t; 9, u), u) x 

(2.11) eipKjds.J J ^ (x(s, 0; 9', J) - x(s, t; 9, u))h(9', u')d9'du' 

respectively. Even if h(9,u>) is not differentiable, we consider Eq. (l2.1ip to be a weak 
solution of Eq. (l2.4p . It is easy in usual way to prove that the integro-ODE ( 12.1(1 has a 
unique solution for any t > 0, and this proves that the continuous model Eq. (12.41) has a 
unique weak solution (12.111) for an arbitrary initial data h(9,u). 

Throughout this paper, we suppose that the initial date /i(6>, u) is of the form 
h(9, u) = h{9)g{u). This assumption corresponds to the assumption for the Kuramoto- 
Daido model (11.31) that initial values {9j(0)}jL 1 and natural frequencies {tOj}jLi are 
independently distributed. This is a physically natural assumption used in many liter- 
atures. In this case, pt(9,u) is written as pt{9,u) = p t (9 \u)g(u), where 

p t (9, u) = h(x(0, t; 9,u)) x 

(2.12) exp[iv- f ds- [ [ %(x(s,0;9',lo') - x(s,t; 9 1 u))h(9')g(u')d9'du' 

L Jo JnJo o9 \ J 

and p t (9,u>) satisfies the same equation as Eq.( l2.4p . 

§ 3. Linear stability of the incoherent solution 



A trivial solution of the continuous model 02.41) . which is independent of 9 and 
t, is given by p t (9,u)) = g(uj)/(2ir), or equivalently /5t(6>,u;) = l/(27r). It is called the 
incoherent solution, which corresponds to the de-synchronized state. Note that in this 
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case r = 0. In this section, we investigate the stability of the incoherent solution and 
the order parameter. 
Let 

/•2tt 

(3.1) Zj(t,u):= e^ je p t {6,u)d6 

Jo 

be the Fourier coefficients of pt(0, to). Then, Zo(t, uS) = 1 and Zj, j = ±1, ±2, • • • satisfy 
the differential equations 

— oo 

The incoherent solution corresponds to the zero solution Zj = for j = ±1, ±2, • • • . 
Since \Zj(t,uj)\ < 1, Zj(t,u) is in the Hilbert space L 2 (R, g(uj)du) for every t : 



\Zj\\h iRig{Lj)dL0) = j \Zj{t,w)\ 2 g{u})dw < 1. 



Thus we linearize the above equation as an evolution equation on L (R, g(ui)doj) 

(3.2) _A = ( j ^i M+j ^T Kfj <p) Zj , j = ±l,±2,--., 

where M. : q(u>) i— > uiq{ui) is the multiplication operator on L 2 (R, g(uj)duj) and V is the 
projection on L 2 (R, g(u))du>) defined to be 

(3.3) Vq{uj) = / q{uj)g{uj)duj. 

Jr 

If we put Po(uj) = 1, V is also written as Vq(uj) = (q,Po), where ( , ) is the inner 
product on L 2 (R, g(u)dw): 



(u)duj. 



(3-4) (91,92):= / 91(^)92(0;)^ 

Note that the order parameter is given as Zf = (Zi, P )- To determine the stability of 
the incoherent solution and the order parameter, we have to investigate the spectrum 
and the semigroup of the operator Tj := j^/— T-M + jy/^lKfjV. 



§ 3.1. Analysis of the operator \f^\M. 

If fj = 0, Tj = j^J—lM. It is known that the multiplication operator M. on 
L 2 (R, g(uj)duj) is self-adjoint and its spectrum is given by a(Ai) = supp(g) C R, where 
su PP(fiO is a support of the density function g. Thus the spectrum of j^—lM. is 

(3.5) a(jV^lM) = jV^l ■ supp(g) = {jV^lX | A G supp(g)} C v^R- 
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The semi-group e J%/ 1 generated by is given as eiv^ Mt q(uj) = e J%/ q(u). 

In particular, we obtain 

(3.6) (e^ Mt qil q 2 ) = [ e^" t q 1 {u)ql&)g{u)du 

Jn 

for any qi , q 2 G £ 2 (R, g(uj)doj). This is the Fourier transform of the function qi(u)q2(u>)g(u>). 
Thus if q± (uj)q2(uj)g{uj) is real analytic on R and has an analytic continuation to a neigh- 
borhood of the real axis, then (ei^ r ~^ Mt qi,q2) decays exponentially as t — > oo, while if 
qi(u>)q2(ui)g(ui) is C r , then it decays as 0(l/t r ) (see Vilenkin [2T]). 
These facts are summarized as follows: 

Proposition 3.1. Suppose that fj = and Eg. \3.2\i is reduced to dZj/dt = 
j\f—\M.Zj. A solution of this equation with an initial value q(u) G £ 2 (R, g(u)du) 
is given by Zj(t) = e^^~^ Mt q{u) = e : ' y/r ~^ ujt q(uj). In particular the linearized order 
parameter Zi(t) = (e^~^ Mt q, Pq) decays exponentially as t — ► oo if g{u) and q{u) have 
analytic continuations to a neighborhood of the real axis. 

The resolvent (A — 1.M) -1 of the operator jy/— 1M. is calculated as 

(3.7) ((X-jV^lM)- 1 qi,q 2 )= [ - ^=^qi{u)^(u)g(u)dtj. 

JnX - jyf-lu) 

We define the function D(X) to be 

(3.8) D(X) = ((A-jv /z T-M)- 1 P , J Po) = / ^^g^duj 

(recall that Pq{oj) = 1). It is holomorphic in C\a(j\/^ 1A4) and will play an important 
role in the later calculation. 

§ 3.2. Analysis of the operator Tj = jy^TM + jy/^lKfjV 

In what follows, we suppose that fj ^ 0. The domain D(Tj) of Tj is given by 
D(.M) fl D(V) = D(A4). Since M. is self-adjoint and since V is bounded, Tj is a closed 
operator [8j. Let g{Tj) be the resolvent set of Tj and cr(Tj) = C\g(Tj) the spectrum. 
Since Tj is closed, there is no residual spectrum. Let cr p (Tj) and cr c (Tj) be the point 
spectrum (the set of eigenvalues) and the continuous spectrum of Tj , respectively. 

Proposition 3.2. (i) Eigenvalues X of Tj are given as roots of 

(3.9) D(X) = . J. , A G C\a(jy/=iM). 

3V-lKfj 
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(ii) The continuous spectrum of Tj is given by 

(3.10) a c {T 3 ) =a{jy/=lM) = jV=l • supp(s). 



Proof. 



(i) Suppose that A G a p (Tj)\a(jy^lM). Then, there exists x G 



L 2 (R, g(uj)duj) such that 

Ax = (jV^TM + jV-LKfjV)x, x^O. 

Since A a(j \/ — 1M.), (A — j\f^lM)~ 1 is defined and the above is rewritten as 

x = (A - jV^lMy^V^lKfjVx 
= jV^lKf J (x,P )(\-jV^lM)- 1 Po(u J ). 

By taking the inner product with Pq(u), we obtain 

(3.11) l=jV=lKf j ((\-jV=lM)- 1 P ,P ) =jV^lKf J D(X). 



This proves that roots of Eq. fl3.9p is in o~ p (Tj). The corresponding eigenvector is given 
by x = (A - j\^lM)~ 1 P (uj) = 1/(A- jV^lu). If A G V-IR, x i L 2 (R, g(u)du). 
Thus there are no eigenvalues on the imaginary axis. 

(ii) This follows from the fact that the essential spectrum is stable under the bounded 
perturbation and that there are no eigenvalues on a(j\/—lji4), see [8]. □ 

§ 3.3. Eigenvalues of the operator Tj and the transition point formula 

Our next task is to calculate roots of Eq. (l3.9p to obtain eigenvalues of Tj = 
jy/^lM + jyf—lKfjV. By putting A = x + \f—\y, x,y G R, Eq. (13.91) is rewritten 

as 



(3.12) 



x 



R x 2 + (ju - y) 2 
ju - y 



g(u>)du 
g(u)du 



Mfj) 

jK\fj\ 2 ' 



In what follows, we suppose that Im(/j) < 0. The case Im(/j) > will be treated in 
Sec. 3. 5. The next lemma is easily obtained. 



Lemma 3.3. 

(i) When Im(/,) < 0, A satisfies Re(A) > for any K > 0. 

(ii) If K > is sufficiently large, there exists at least one eigenvalue A near infinity. 
(Hi) If K > is sufficiently small, there are no eigenvalues. 
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Proof. Part (i) of the lemma immediately follows from the first equation of Eq. ( 13.121) . 
To prove part (ii) of the lemma, note that if |A| is large, Eq. ( 13.91) is rewritten as 

1 + o( 1 ' 



Thus the Rouche's theorem proves that Eq. (I3.9H has a root A ~ j^—lKfj if K > 
is sufficiently large. To prove part (iii) of the lemma, we see that the left hand side of 
the first equation of Eq. ( 13.121 ) is bounded for any x,y G R. To do so, let G(u) be the 
primitive function of g(oj) and fix 5 > small. The left hand side of the first equation 
of Eq. (13. 1 2[) is calculated as 

xg(u)du 



R x 2 + (ju-y) 2 



xg(u)du | f y '^ 5 xg(u)du> ^ ^v/j+ 5 xg(u)du 



in 



/ j+s x 2 + (juj - y) 2 i-oo x 2 + (ju - y) 2 J y /j-s x 2 + (ju - y) 2 
xg(u)du ^ j-y/j-s xg(u)du 



\2 



y/j+s x 2 + (ju - y) 2 J_ 00 x 2 + (ju - y) 

(Gfe/i + *) " G(v/j - *)) + J+t-tw ^ 

The first three terms in the right hand side above are bounded for any x, y G R. Since 
G is continuous, there exists a number £ such that the last term is rewritten as 

This is bounded for any x, y G R. Now we have proved that the left hand side of the 
first equation of Eq. (13.12p is bounded for any x > 0, although the right hand side tends 
to infinity as K — > +0. Thus Eq. ( 13.91) has no roots if K is small. □ 

Lemma 3.3 shows that if K > is sufficiently large, the trivial solution Zj = 

of the system dZj/dt = TjZj is unstable because of the eigenvalues with positive real 

(i) 

parts. Our purpose in this subsection is to determine the bifurcation point K c , which 

(i) 

is the minimum value of K such that if K < Kc , the operator Tj has no eigenvalues on 
the right half plane. To calculate eigenvalues A = \(K) explicitly is difficult in general. 
However, note that since zeros of a holomorphic function do not vanish because of the 
argument principle, \(K) disappears if and only if it is absorbed into the continuous 
spectrum a(j\/—lA4), on which D(X) is not holomorphic. This fact suggests that to 
determine it is sufficient to investigate Eq. (l3.9p or Eq.( )3.12]) near the imaginary 

axis. Since we are interested in X(K) absorbed into a(j^/—lA4) C \/— 1R, take the 
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supp(g) 

Figure 3. A schematic view of behavior of roots A of Eq. ( 13.91) when K decreases. 



limit x -»• +0 in Eq. (l3TT2]) : 



(3.13) 



x 



lim , ., 

x ^+°Jrx + KJU - y) 2 

< x ^+°JrX 2 + (juj - y)' 



g(u)du 
g(u)dw 



Re(/ 3 Q 



These equations determine K n and y n such that one of the eigenvalues A n (X) converges 
to \/—ly n as -ftT — > ii^ + (see Figj3]). To calculate them, we need the next lemma. 



Lemma 3.4. (i) Suppose that X n (K) 
continuous at u = y n . 
(ii) If g(u>) is continuous at u> = y, then 



-ly n as K — > K n . Then, g(u) is 



(3.14) 



lim 



x 



+o J R x 2 + (ju -y) 2 



g(u)du = ng(y/j)/j. 



Proof. To prove (i), suppose that g(u) is discontinuous at u> = without loss of 
generality. 

STEP 1: At first, we suppose that g(u) is piecewise continuous. Put g(+0) = 
h+, g(—0) = h- and h + ^ h_. In this case, for any e > 0, there exists 5 > such 
that if — 5 < u < 0, then \g(u) — h-\ < e and if < u < 5, then \g(u) — h+\ < e. For 
Eq. ( 13.91) . we suppose |A| = \x + \^Ty\ < 5 and y > 0. The case y < is treated in a 
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similar manner. We calculate D(X) as 



D(X) 



5 A - jy/-lu 



-dw + 



-doo 



+ 



1 



o X-jV-luJ 



(g(oo) — h + + h + ) doo + 



x-r 



-loo 



(g(oo) — h- + h-)dw 



+ 



x-r 

s 



-loo 



doo + 



X+ jy/-l0O 



du 



1 



o A-j 

s 



-loo 



(g(oo) — h + ) doo + 



1 



(3.15) +/ 1 (h+ - h-) duj + h- [ 

Jo X-jV-loo Jo 



o A + jV-lw 
doo 



(g(—u) — h-)duo 



x-r 



-loo 



+ h. 



doo 



o X + jV-loo 



Since |A| < 5, there exists a positive number M, which is independent of A, such that 

1 g(±v) 



Thus |-D(A)| is estimated as 

\D(X)\ > \h+ — h-\ [ 
Jo 



s A=fjV-1w 



doo 



-doo < M. 



h- 
j 



log 



A - jV^luj 



doo 



-IX- jS 



'o |A — 
- 2M. 



■lu\ 



doo 



o \\ + jV-lv\ 



Since y > 0, |A — jy/— loo\ < \X + jy/— Too\. This shows that 



(3.16) \D(X)\ > (\h+ -h-\ -2e) 



doo 



o \X-jV-loo\ 



h- 

— log 

J 



-1A - jd 



2M. 



The right hand side tends to infinity as A — > if 2e < \h + — h—\. This proves that 
Eq. ( 13.91) has no roots at A = for positive K. 

STEP 2: In general, since g(oo) is a non- negative measurable function, there exists 
a monotonic increasing sequence {g n (oo)}^ =1 of non- negative simple functions such that 
g n (u>) — > g(oo) for each oo. In particular if g(oo) is discontinuous at u> = 0, we can choose 
{g n (u)}^ =1 so that g n (oo) is discontinuous at oo = for any n G N. Then, the proof is 
done in the same way as STEP 1 by approximating g(oo) by g n (u). 

(ii) The formula Eq. ( 13.141) is proved in Ahlfors pQ. □ 

Let (y, K) be one of the solutions of Eq. (13. 13ft . Since g(oo) is continuous at oo = y, 
substituting it into the first equation of Eq. ( 13.130 yields 



(3.17) 



ng(y/j) 



K\f: 



■ |2 ' 
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Substituting K = —Im(fj)/(Tr\fj\ 2 g(y/j)) obtained from the above into the second 
equation of Eq. ( 13.131) results in 

(3-18) lim / 2 ^.~ y . 2 g(u)du = ^Mg(y/j). 

This equation for y determines imaginary parts to which X(K) converges as Re(A(_ftT)) — > 
+0. Let yi,y 2 , ■ ■ ■ be roots of Eq. fl3TT8jl . Then, 

(3-19) K n = |7| 2 m i /j) /-V n = l ^-'- 

give the values such that Re(A(_KT)) — > as K —>■ K n + 0. Now we obtain the next 
theorem. 



Theorem 3.5. Suppose that Im(fj) < 0. Let y\,yi-,--- be roots of Eq.fiSHE)- 

Put 

(3.20) K^i=w£K r ,- ~~ lm(fj) 



L n 

n 



Mfj\ 2s ^Pn9(yn/j)' 



(i) 

If < K < Kc , the operator Tj has no eigenvalues, while if K is slightly larger than 
(i) 

K c , Tj has eigenvalues on the right half plane. 

Note that inf n K n is positive because of Lemma.3.3 (iii). As a corollary, we obtain 
the transition point (bifurcation point to the partially locked state) conjectured by Ku- 
ramoto |10|: 



Corollary 3.6 (Kuramoto's transition point). Suppose that the probability den- 
sity function g(u) is even and maxg(u) = g(0). If Re(fi) = and Im(/i) = —1/2 (it 
corresponds to f(9) = sin 9 in Eq. U.3\) ). then K^ defined as above is given by 



(3-21) 




When K > Kc , the solution Z\ = of Eq. \3. 2\) is unstable. 

§3.4. Semi-group generated by the operator Ti (Im(/i) < 0) 

Since we are interested in the dynamics of the order parameter Z®{t) = (Zi,P ), 
in what follows, we consider only j = 1 while cases j = 2, 3, • ■ ■ are investigated in 
the same way. Theorem 3.5 shows that Kc is the least bifurcation point and the 
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trivial solution Zi(t,u) = of Eg. (I3.2p is unstable if K is slightly larger than If 
< K < the spectrum of T\ is on the imaginary axis: cx(Ti) = a(yf—\M), and 

thus the dynamics of Z\ is nontrivial. In this subsection, we investigate the dynamics 
of Z\ and the order parameter for < K < We will see that the order parameter 

may decay exponentially even if the spectrum lies on the imaginary axis because of 
existence of resonance poles. 

Since \f-\M. has the semi-group e ^~^ Mt and since V is bounded, the operator 
T\ = y/—lM. + \/—lKf\P also generates the semi-group (Kato [8]), say e Tlt . A solution 
of Eq. (13.2p with an initial value q(u>) G L 2 (R, g(u)duj) is given by e Tlt q(u). The e Tlt is 
calculated by using the Laplace inversion formula 



(3.22) 



lim 

y^oo 27Ti 



1 



'111 



e A *(A-T 1 )" 1 rfA, 



where x > is chosen so that the contour is to the right of the spectrum of Ti 
(Yosida [22]). At first, let us calculate the resolvent (A — Ti) -1 . 



Lemma 3.7. For any q(uS) G L 2 (R, g(uj)dui), the equality 



(3.23) 
holds. 



1 - ^lKhD(X) 



Proof. Put 

R(X)q := (A - Ti)" 1 ? = (A - V^IM - V^lKhV^q, 

which yields 

(A - V^lM)R(X)q = q + V^lKf 1 VR(X)q 

= q + y^TKf 1 (R(X)q,P )Po. 

This is rearranged as 

R(X)q = (A - ^\M)- l q + V=lKf 1 {R(X)q, P )(A - sTIM^Pq. 
By taking the inner product with Pq, we obtain 

(R(X)q, P ) = ((A - V^lMy'q, P ) + y/=lKf 1 (R{X)q, P )D(X). 
This proves Eq. (ET25j) . 



□ 
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Let Z®(t) = (Zi,Pq) be the order parameter with the initial condition Zi(0,u) = 
q(u). Eqs.(H£22D and (13T23D show that Z$(t) is given by 

(3,4, m = (e- ? ,Po) = H m r^> ((A-W'^) dA 

One of the effective way to calculate the integral above is to use the residue theorem. 
Recall that the resolvent (A — Ti) _1 is holomorphic on C\a(T±). Since we assume that 
< K < K { c\ Ti has no eigenvalues and the continuous spectrum lies on the imaginary 
axis : a{T x ) = a(y/=lM) = ■ supp(#). Thus the integrand e Ai F (A) in Eq.(E21D 
is holomorphic on the right half plane and may not be holomorphic on er(Ti). However, 
under assumptions below, we can show that -Fb(A) has an analytic continuation F±(X) 
through the line er(Ti) from right to left. Then, i*i(A) may have poles on the left half 
plane (the second Riemann sheet of the resolvent), which are called resonance poles |17j . 
The resonance pole \x affects the integral in Eq. (l3.24p through the residue theorem (see 
FigJ3D • In this manner, the order parameter Z®(t) can decay with the exponential rate 
Re(^u). Such an exponential decay caused by resonance poles is well known in the theory 
of Schrodinger operators [17] , and for the Kuramoto model, it is investigated numerically 
by Strogatz et al. [19] and Balmforth et al. [2]. 

At first, we construct an analytic continuation of the function Fq(\). 



Lemma 3.8. Suppose that the probability density function g(uj) and an initial 
condition q{u) are real analytic on R. If g(uj) andq{ui) have meromorphic continuations 
g*(X) and q*(X) to the upper half plane, respectively, then the function -Fo(A) defined 
on the right half plane has the meromorphic continuation Fi(X) to the left half plane, 
which is given by 

((A - y^lMy'q, P ) + 2nq*(-^T\)g*(-V^l\) 



(3.25) 



Fi(A) 



1 



-lKf 1 D(X)-27r^Kf 1 g*(- 



-1A) 



Proof. By the formula ( 13. 14j) . we obtain 
(3.26) 

lim ((A-v^-MrV^o)- hm ((X-V^lM)~ 1 q,Po)=2nq(Im(X))-g(Im(X)). 

Re(A)^+0 Re(A)^-0 

Thus the meromorphic continuation of ((A — V~ lM)~ l q, Pq) from right to left is given 

by 



((A 
((A 



(3.27) 

This proves Eq. f[3T25j) . 



-lMy^Po) 
-lM)- 1 q,P ) + 2nq*(- 



-lX)g*(- 



(Re(A) > 0), 
TA) (Re(A) < 0). 



□ 



16 



Hayato Chiba 



Poles of Fi(X) (resonance poles) on the left half plane are given as roots of the 
equation 

1 



(3.28) 



D(X) + 2ng* 



-1A) 



Re(A) < 



and poles of the function q*( — \J — 1A). In the next theorem, we suppose for simplicity 
that q*(—\/—lX) has no poles. Now we calculate the order parameter Z^it). 



Theorem 3.9. For Eg. ( Iff. 2\) with j = 1, suppose that 

(i) Im(/i) < and < K < Kp . 

(ii) the probability density function g(iv) is real analytic on R and has a meromorphic 
continuation g*(\) to the upper half plane. 

(Hi) an initial condition q{u) is real analytic on R and has an analytic continuation 
g*(A) to the upper half plane. 

(iv) there exists a positive number 5 such that |-Fi(A)| — > as |A| — > oo in the angular 
domains 

(3.29) |arg(A)| < 5, |arg(A) - ir\ < 6. 

(v) there exist positive constants D and f3 such that 

(3.30) |Fi(A)| <De m 
in the angular domain tv/2 + 5 < arg(A) < 3tv/2 — d. 

Then, there exist resonance poles ofTi on the left half plane. Let cti, «2, ■ ■ ■ be resonance 
poles such that \a\\ < \a%\ < • ■ ■ . Then, there exists a positive constant t Q such that the 
order parameter is given by 

oo 

(3.31) Z°(t) = (e T ^q,P ) = J2Pn(t)e ant , t > t 

n=l 

where p n (t) is a polynomial in t. In particular, Z®(t) decays exponentially as t — >• oo. 

Proof. At first, we prove the existence of resonance poles. Resonance poles are 
roots of Eq. ( 13.28] ). which is the analytic continuation of the equation (13.9P for j = 1. 
Thus one of the resonance poles is obtained as a continuation of an eigenvalue X(K). 
Recall that X(K) converges into the imaginary axis as K — > + 0. To prove that 
there exists a resonance pole on the left half plane when K < Kc , we have to show 
that X(K) does not stay on the imaginary axis for K < KP. Differentiating Eq. (l3T9j) 
with respect to K, we obtain 



(3.32) 
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Figure 4. The contour for the Laplace inversion formula. 



which proves that X'(K) 7^ 0. Further, roots y of Eq. (l3.18|) . which determines eigenval- 
ues on the imaginary axis, are isolated because both side of Eq. ( 13.181 ) are analytic with 
respect to y. This means that X(K) can not move along the imaginary axis. This proves 
that an eigenvalue X(K) gets across the imaginary axis from right to left as K decreases 
from Kc , which gives a root of Eq. (13.28p . Note that there may exist resonance poles 
which are not continuations of eigenvalues (see Example 3.11). 

Next, let us prove Eq. ( 13.311) . Let d > be a small number and r sufficiently large 
number. Take paths C\ to Cq as are shown in FigJU 



C\ = {d + V—ly I — r < y < r}, 

C 2 = {x + V^lr \0<x<d}, 

C 3 = {re^ ld \n/2<9< n/2 + 5}, 

C 4 = {re^ 16 \n/2 + S<6< 3tt/2 - 5}, 



and C5 and Cq are defined in a similar way to C3 and C2, respectively. We put C = 

Let «i, «2, • • • , ct n (C) be resonance poles inside the closed curve C, where we assume 
that there are no resonance poles on the curve C by deforming it slightly if necessary. 
Let Ri(t), Rzit), ■ ■ ■ ,-R n ( c )(t) be corresponding residues of e xt Fi(X), respectively. Note 
that if ctj is a pole of -Fi(A) of order rrij, Rj(t) is of the form Rj(t) = pj{t)e a:it with a 
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polynomial Pj(t) of degree rrij — 1. By the residue theorem, we have 

n(C) 

2nV^T^ Rj(t) = e xt F (X)dX + e xt F 1 (X)d\. 

j=i JC6+C1+C2 J C3+C4+C5 

The integral J Ci e At Fo(A)<iA/(27r\/^T) converges to as r — > 00. It is easy to show 

that the integrals along C2, C3, C5, C*e tend to zero as r — >■ 00 because of the assumption 
(iv). We have to estimate the integral along C4 as 



e Fi(X)dX 

c 4 



< / re rtcos0 ^{re^ 16 )^ 

J-k/2+5 

< max \F 1 (re^ Td )\ / re rtcosd d6 

7r/2+5<e<S7r/2-S Jn/2+5 
■ it/2 

<De pr / 2re- rtsin ^d<b 



r/2 



(3.33) < De? r ■ j {e-^h 

Thus if t > to := max {,3, n(3/(25)}, this integral tends to zero as r — > 00. □ 

Example 3.10. If g(u;) is a rational function, the assumptions are satisfied when 
q* (A) is bounded on the upper half plane. In this case, the number of resonance poles 
is finite and thus Eq. (13.3ip becomes finite sum. For example if g(u) = l/(7r(l + u 2 )) is 
the Lorentzian distribution, a resonance pole is given by A = \/ —\Kf\ — 1 (a root of 
Eq. ( 13.28( 1). Therefore Z\{€) decays with the exponential rates Re(A/ — l-^/i — 1)- 

Example 3.11. If g(oS) is the Gaussian distribution, the assumptions are sat- 
isfied when q*(X) is of exponential type; that is, there exist positive constants C and 
(3 such that |?*(A)| < Ce^L Since the analytic continuation g*(X) has an essential 
singularity at infinity, there exist infinitely many resonance poles and they accumulate 
at infinity. 



§3.5. Semi-group generated by the operator T\ (Im(/i) > 0) 

In Sec. 3.1 and Sec. 3. 4, we investigate the semi-group generated by the operator 
T\ = \J-\M. + \J—\Kf\P for the cases f\ = and Im(/i) < 0, respectively. In this 
subsection, we consider the case Im(/i) > 0. 

Theorem 3.12. Suppose that the assumptions (ii) to (v) of Thm.3.9 hold. If Im(/i) > 
0, for an arbitrarily fixed K > 0, the order parameter Z®{t) = (e Tlt q,Po) decays 
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exponentially as t — ► oo. 

We show an idea of the proof. If Im(/i) = 0, T\j\f—\ = M. + KRe(fi)V is self-adjoint 
and a rank one perturbation of the multiplication M. . By Theorem X-4.3 in [8], T\j\J—\ 
and M. are unitarily equivalent. Since {e/~^ Mt q, Pq) decays exponentially (see Sec.3.1), 
we can prove that so is (e Tlt q, Pq). 

If Im(/i) > 0, change the parameter as K i— ► — K. Then, the problem is reduced 
to the case K < and Im(/i) < 0, and Thm.3.12 is proved in a similar manner to the 
proof of Thm.3.9. 
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